The gravitational redshift of boson stars 
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We investigate the possible gravitational redshift values for boson stars with a self-interaction, 
studying a wide range of possible masses. We find a limiting value of zn m ~ 0.687 for stable 
boson star configurations. We compare theoretical expectation with the observational capabilities 
in several different wavebands, concluding that direct observation of boson stars by this means will 
be extremely challenging. X-ray spectroscopy is perhaps the most interesting possibility. 
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I. INTRODUCTION 

The idea of the boson star goes back to Kaup in 1968 
HJ. A boson star is a gravitationally bound collection 
of bosonic particles, arising as a solution of the Klein- 
Gordon equation coupled to general relativity. Many in- 
vestigations of the possible configurations have been car- 
ried out; for reviews see Ref. For non-self-interacting 
bosons of mass m, the mass of a typical configuration is of 
order mpj/m, to be compared with a typical neutron star 
mass of Tipj/m^g^Qjj, which is about a solar mass. Here 
topi is the Planck mass. Unless the bosons are extremely 
light, boson stars made from non-interacting bosons are 
orders of magnitude lighter than typical stars. 

The situation is very different if the boson stars have 
even a very weak self-interaction; for a first investiga- 
tion of self-interactions see Ref. Colpi, Shapiro and 
Wasserman Q showed that the maximum mass of stable 
configurations is then of order A^TOpj/rra 2 , where A is 
the scalar field self-coupling, normally assumed to be of 
order unity. Then boson star configurations exist with 
mass (and radius) similar to that of neutron stars, if the 
bosons, like neutrons, have a mass around 1 GeV. They 
can also be much heavier, should the bosons be lighter. 
In this paper, we allow ourselves to consider a very wide 
range of possibilities for the boson star mass and radius. 

It remains unknown whether or not boson stars could 
actually form in our own Universe, for example through 
a gravitational instability process, though there has been 
some discussion in the literature |J . The most favourable 
situation would be for the dark matter in the Universe to 
be in the form of very weakly interacting bosons, a sub- 
stantial fraction of which manage to form boson stars; 
such a process may be aided if the bosons have a rea- 
sonably strong coupling to each other, but not to con- 
ventional matter. If boson stars can form, they provide 
an alternative explanation for stellar systems in which 
an object is inferred to have a high mass; convention- 
ally, a 'star' with mass greater than a few solar masses is 



assumed to be a black hole (since this is above the max- 
imum permitted mass for a neutron star), but a boson 
star offers a more speculative alternative. 

Dark matter is characterized by its extremely weak in- 
teraction with conventional matter; even though 90% of 
the mass of the galaxy is in this form, it is not directly 
visible. In this paper, we investigate the implications of 
assuming that the material from which boson stars are 
made is similarly extremely weakly interacting, so that its 
only interaction with neighbouring baryonic material and 
photons is gravitational, just as the relation between the 
visible galaxy and its dark matter halo.0 However, while 
a galaxy halo can be described using Newtonian theory, 
boson stars close to the maximum allowed mass are gen- 
eral relativistic objects. This gives a new characteristic 
of such objects, a gravitational redshift. In this paper, 
we examine whether radiating baryonic matter moving in 
the gravitational potential of a boson star could be used 
as an observational signal of boson stars. 



II. THE EINSTEIN-KLEIN-GORDON 
EQUATIONS 

The Lagrange density of a massive complex self- 
gravitating scalar field reads 



C 



'^R 

8tt 



C/(|$| 2 ) 



(1) 



where R is the curvature scalar, g the determinant of 
the metric g^ v , and $ is a complex scalar field with a 
potential U. We take ft = c = 1. Then we find the 
coupled system of equations 



*An example already existing in the literature is the boson- 
fermion star which is made up of bosons and neutrons 
interacting only gravitationally. 
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is the energy-momentum tensor and 
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the generally covariant d'Alembertian. Within the model 
we want to have an additional global f (1) symmetry, so 
we can take the following potential 
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where m is the scalar mass and A a dimensionless con- 
stant measuring the self-interaction strength. 

The method for finding solutions is well known Elpl. 
For spherically symmetric solutions we use the static line 
element 



ds 1 



e u{r) dt 2 - e^ r) dr 2 - r 2 (dd 2 + sin 2 6 dip 2 ) . (7) 



The most general scalar field ansatz consistent with this 
metric is 



$(r,t) = P(r)e~ 



(8) 



where u> is the frequency. Notice that to find finite mass 
solutions, the scalar field must carry a time-dependence 
of this form, which leaves the energy-momentum tensor, 
and hence the metric, time-independent. 

The non-vanishing components of the energy- 
momentum tensor are 

T Q ° = p= X -\u?P 2 {r)er v + P l2 (r)er» + U] , (9) 

Ti 1 = Pr = l -[u 2 P 2 {r)e-» + P' 2 (r)e-" - U] , (10) 
T 2 2 = T 3 3 = p ± with 

P± = -^[u 2 P 2 (r)e- v - P' 2 {r)e-» - U] (11) 

where ' = dj dr. Note that the pressure is anisotropic; 
there are two equations of state p r = p — U and p± = 
p-U - P' 2 (r)e-v. 

The decisive non- vanishing components of the Einstein 
equation are 
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the two further components giving degenerate equations 
because of the Bianchi identities. 

The differential equation for the scalar field is 
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For the rest of our paper, we employ the dimension- 
less quantities x — mr, a — y/An P/mpi and A = 
ATOpj/47rm 2 . In order to obtain solutions which are regu- 
lar at the origin, we must impose the boundary conditions 
ct'(0) = and p{0) = 0. 

Once the fundamental parameter A has been fixed, the 
only remaining freedom is the value of the scalar field 
at the centre of the star (the scalar field mass simply 
contributing an overall scaling); the initial value ^(0) is 
determined by the boundary condition of Minkowskian 
spacetime at infinity. It is well known that as the cen- 
tral scalar field value is increased, initially the inferred 
mass of the star (determined from the asymptotic form 
of the metric) increases up to a maximum, after which 
it exhibits oscillatory behaviour, shown in Figure |l]. The 
configurations up to this maximum mass are dynami- 
cally stable, but those beyond it are unstable [f|-|U|] and 
should not be considered physical. 



III. GRAVITATIONAL REDSHIFT 

In this Section we calculate the redshift produced by 
the gravitational potential of a boson star. Under our 
assumption that the scalar particles have no interaction 
with the baryonic matter other than gravitationally, the 
baryonic matter can penetrate up to the center of the 
boson star. If the baryonic matter emits or absorbs ra- 
diation within the gravitational potential, the spectral 
feature will be redshiftcd. 

The gravitational redshift z within a static line element 
is given by (see e.g. Ref. Jl0||) 



1 + z„ = 



(15) 



where the emitter and receiver are located at P; n t and 
i? ox t respectively. The maximum possible redshift for a 
given configuration is obtained if the emitter is exactly at 
the center i?i nt = 0. The receiver is always practically at 
infinity, hence exp [^(i? cxt )] = 1. For all other redshifts 
in between, we define the redshift function 



1 + z g (x) = exp 



u(x) 



(16) 



Table [l| gives the results for the maximum mass boson 
star for different A. The maximum mass gives the highest 
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value one can obtain from stable configurations (unsta- 
ble configurations can yield very high redshift values — 
we show one as an example). One recognizes that with 
increasing self-interaction values A also the maximal red- 
shift value grows. For unstable states very high z g values 
can be found. For A — > oo, we find the asymptotic value 
z\ im ~ 0.687, as shown in Fig. ||. 

Figure || shows the decrease of the asymptotically mea- 
sured redshift as a function of radius, for different A val- 
ues. 

The highest central value is quite interesting. In the 
early years of quasar detection, it was speculated whether 
the origin of the redshift could be gravitational rather 
than cosmological; see, for example, Refs. |lj|-|lj|. The 
Schwarzschild interior solution for a perfect fluid sphere 
with constant density gives a maximal redshift on the 
surface of 2 (cf. 14 1), but this matter form is unrealistic 
because it has a transmission velocity exceeding c. The 
value 2 is also an upper limit for all static fluid spheres 
whose density does not increase outwards @. Bondi Jl|] 
applied constraints on a perfect fluid with isotropic pres- 
sure and looked for stable objects; from these results, a 
maximal surface redshift of z m ax — 0.62 can be derived, 
a value which is quite close to our central one of 0.687 
for the anisotropic boson matter fluid. For a neutron 
star, one finds a central value of 0.47 p6| , whereas for a 
star consisting of two kinds of fermions one has a max- 
imal central redshift (for a stable configuration) of 1.22 
|L7| . Different forms of matter result in completely differ- 
ent central redshift values. Just a small self- interaction 
A = 1 for the bosonic particles gives a redshift value 
about 0.49, just above that of the neutron star. 

In general, observed redshift values could consist of a 
combination of cosmological and gravitational redshift. 
Their relation is given by 



l + z= (1 + z c )(l + z g ) 



(17) 



because of the relativistic velocity addition formula. The 
objects with highest redshift values today are the quasars 
with redshifts up to about five. If one assumes that they 
would have an incorporated gravitational redshift, and 
our solutions give a maximal value of about 0.7, then we 
have a cosmological redshift of 2.53. Even if for some 
reason they possessed a maximal gravitational redshift, 
this could not entirely explain the observed redshifts. In 
any case, there are no indications of a non-cosmological 
redshift for quasars at present. 



IV. DETECTABILITY OF BOSON STARS 

The mass M of a boson star composed of non-self- 
interacting particles is inversely proportional to m, while 
the mass of a self-interacting boson star is proportional to 
V~\/m 2 M. Taking A ~ 1, then for small m (to be pre- 
cise, provided A> 1) the slef-interacting star is much 



more massive. For example, if we want to get a boson 
star with a mass of order 10 33 g (a solar mass), then we 
need m ~ 10~ 10 eV for A = 0, or m cx A 1 / 4 GeV if 
A 3> 10~ 38 (we see that the self-coupling has to be ex- 
traordinarily tiny to be negligible). In this example, the 
scalar particle has a mass comparable to a neutron, lead- 
ing to a boson star with the dimensions of a neutron star. 
If we reduce the scalar mass further, to m ~ 1 MeV, then 
we find M ~ 10 39 \/A g and R ~ 10 6 VA km; this radius 
is comparable to that of the sun, but encloses 10 6 solar 
masses. These parameters are reminiscent of supermas- 
sive black holes, for example as in Active Galactic Nu- 
clei; the mass-radius relation is effectively fixed just by 
the objects being relativistic. An exceptionally extreme 
example is to take m ~ 1 eV, giving M ~ 10 18 -\/AMq 
and R ~ lO 3 ^ parsecs. 

These numbers must be compared with the critical 
matter density in the Universe, around 10 12 M Q Mpc -3 . 
This sets the scale of the typical distance expected to the 
nearest object of a given mass, provided such objects are 
common. The largest gravitationally bound objects in 
the present Universe are rich galaxy clusters, with masses 
of order 10 15 Mq, which require the assembly of material 
from a volume around 10 Mpc on a side. 

In all cases, the density of the boson stars makes their 
direct detection difficult; in particular, they cannot be 
resolved in any waveband. The best optical spectroscopy 
has a resolution of order of an arcsecond, nowhere near 
good enough. The best imaging resolution is in the radio 
band; for example the VSOP program p8| , which in- 
cludes the recently launched Japanese HALCA satellite, 
offers resolution of around a milli-arcsecond, but this ap- 
proach is limited because spectroscopy cannot be carried 
out. At a distance of lOMpc, we could resolve a boson 
star with a radius of O.lpc, corresponding to a mass of 
10 15 M Q , but rich clusters cannot possibly be dominated 
by a single massive object of that sort. Any reduction of 
the boson star mass requires it to be closer, in propor- 
tion to the mass M, but the amount of mass available 
in that volume is decreasing as the cube of the distance. 
The situation therefore rapidly becomes worse at closer 
distances; the nearest object of a given size is expected 
to be far too distant to be resolvable. For example, to be 
resolvable at a distance of one kiloparsec, a boson star 
would need a radius of 10 13 cm (roughly the Earth-Sun 
separation) requiring a mass of 10 11 Mq, almost the mass 
of an entire galaxy. 

However, even if boson stars cannot be directly re- 
solved, their influence might still be visible if material 
in their vicinity is sufficiently luminous. It is necessary 
to find a certain amount of luminous matter within the 
gravitational potential of the boson star. This could, 
for example, be HI gas clouds as seen in galaxies. From 
the rotation curves calculated below, one imagines that 
the hydrogen gas would probably be excited, hence HII, 
because of the kinetic energy gain in the gravitational 
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potential of the boson star. One might also expect accre- 
tion discs about boson stars, though there the luminosity 
could be dominated by regions outside the gravitational 
potential and the boson star would be indistinguishable 
from any other type of compact object. One should men- 
tion that bright HII regions are found within the host 
galaxies of quasars p~S| ] . 

A more dramatic possibility would be to house a com- 
plete nuclear burning star within the gravitational po- 
tential of the boson star. Unfortunately that does not 
seem feasible, since we have already seen that to have 
a solar radius the boson star must weigh 10 6 M Q . A 
more generic situation therefore would be for lighter bo- 
son stars to be contained completely within conventional 
stars, influencing the stellar structure by means of their 
additional gravitational interaction. Such a situation has 
in fact been studied for a degenerate neutron core inside 
a supergiant, the so-called Thorne-Zytkow objects [ po[ . 
If the bosons are completely non-interacting, then the 
boson star might indeed be able to survive such an en- 
vironment. With an interaction present, one would have 
to examine whether or not the boson star could be dis- 
sociated by high energy photons in the stellar interior, 
so that instead the scalar particles circulate freely within 
the star. Historically, such a situation was proposed as a 
possible resolution of the solar neutrino problem, under 
the name cosmion 

The results from the last section showed that a boson 
star, like a conventional star, cannot be resolved in it- 
self but rather would have to be part of a larger area. 
We know that the luminosity from the boson star will 
be redshifted. Hence we can distinguish between the lu- 
minosity of the boson star and other areas if there is a 
particularly strong spectral line. The luminosity Lboson 
could be derived from the redshifted tail of the emission 
line subtracting the 'normal' part. First, we assume that 
the gas distribution is everywhere homogeneous within 
the resolved area. This could be, in first order, valid for 
almost Newtonian boson stars. Let L be the complete 
measured luminosity; then we have L — L noima \ + L-^ oson , 
where we divided L into its two parts of the 'normal' non- 
redshifted luminosity and the one from the boson star 
area. Both parts are proportional to the product of the 
relation A of the boson star area to the total resolved 
area and the total luminosity L, hence, Lboson = ^-L 
while Lnormai = (1 — A)L. The calculated value of A 
gives only an upper limit for the size of the boson star. 

Given the inability to resolve boson stars, the most 
promising technique is to consider a wave-band where 
they might be extremely luminous. The best example is 
in X-rays. A very massive boson star, say 1O 6 M is likely 
to form an accretion disk and since its exterior solution is 
Schwarzschild it is likely to look very similar to an AGN 
with a black hole at the center. In X-rays it may be pos- 
sible to probe very close to the Schwarzschild radius; it 
has been claimed by Iwasawa et al. p2] that using ASCA 



data they have probed to within 1.5 Schwarzschild radii. 
Because this is inside the static limit for a non-rotating 
black hole, they conclude that a Kerr geometry is re- 
quired. We note that a boson star configuration provides 
a more speculative alternative, giving a non-singular so- 
lution where emission can occur from arbitrarily close to 
the center. The signature they use is a redshifted wing of 
the Iron K-line, indicative of emission from deep within 
the gravitational potential well. If such techniques have 
their validity confirmed, it may ultimately be possible 
to use X-ray spectroscopy to map out the shape of the 
gravitational potential close to the Schwarzschild radius 
or boson star, and perhaps to differentiate between the 
theoretical predictions for a Kerr black hole and for a 
rotating boson star. 

V. ROTATION CURVES 

We end by considering the rotation curves about a 
boson star in more detail. For the static spherically 
symmetric metric considered here circular orbit geodesies 
obey 

vl = rv'e v /2 = e v (e» - 1)12 + ^- Pr r 2 e^ +U /2 

m P1 

~ M(r)/r + -^lp r rV + 72 , (18) 

which reduces for a weak gravitational field into the New- 
tonian form v N cwt = M(r)/r if p r = 0. 

Figure || shows the rotation curves for the cases A = 
and 300. The curves increase up to a maximum followed 
by a Keplerian decrease. The possible rotation velocities 
circulating within this gravitational potential are quite 
remarkable. Their maximum reaches more then one-third 
of the velocity of light; if boson stars can form, then such 
enormous rotation velocities are not necessarily signa- 
tures of black holes. The matter possesses an impressive 
kinetic energy of about 6% of the rest mass. If we were 
to suppose that each year a mass of 1M Q transfers this 
amount of kinetic energy into radiation, a boson star had 
a luminosity of 10 44 erg/s. 

Figure |5| shows how the rotation velocities decrease 
with decreasing initial value <r(0) at fixed A. Table || 
gives the maximal velocity. For Newtonian solutions the 
rotation velocity is low and quite constant over a larger 
interval. 

VI. DISCUSSION 

In studying a wide range of possible masses and cou- 
pling for the bosons, we have unveiled a range of possi- 
bilities for boson stars to have an observational impact. 
In the literature it is common to consider boson stars to 
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be like neutron stars, but this relies on specific choices 
for the boson mass for which there is no strong motiva- 
tion. A boson star could have quite different dimensions 
both in size and mass. In analogy to the dark matter 
halos of galaxies, a boson star could also be transparent, 
i.e. without any electro- weak and strong interactions. Be- 
cause of the general relativistic background of the boson 
star solutions, a gravitational redshift has to be taken 
into account. We find that for stable boson stars a cer- 
tain redshift value cannot be exceeded. This is in good 
agreement with results for general fluids and especially 
the neutron star model where one finds similar redshift 
values. We mention that if one considers instead bosonic 
particles without mass or self-interaction, one finds solu- 
tions which can be applied to fit rotation curve data of 
spiral and dwarf galaxies, and which may also give very 
high redshift values fl23| . 

It will be observationally very challenging to use the 
gravitational potential of these configurations to provide 
an observational test. Although radio offers by far the 
best imaging resolution, it seems that it may be spec- 
troscopy using X-rays which offers the best opportu- 
nity to probe the central regions of strong gravitational 
sources. 
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TABLE I. The redshift values 

z for different self-interaction values A and initial values of 
the scalar field a. The addition (max) means the maximal 
central density up to which stable boson star states exist for 
the choice of A. 



_A ajfi) z 

^20 0.067 0.0656984 (max) 

-15 0.086 0.0869732 (max) 

-10 0.122 0.129636 (max) 

-5 0.195 0.237231 (max) 

0.025 0.0315 

0.271 0.456516 (max) 

1.95 112.14 

5 0.253 0.565177 (max) 

10 0.225 0.610195 (max) 

20 0.184 0.642457 (max) 

30 0.158 0.652836 (max) 

50 0.128 0.665739 (max) 

75 0.107 0.673652 (max) 

100 0.0935 0.674607 (max) 

200 0.0673 0.681385 (max) 

300 0.0552 0.682227 (max) 



A > 1 log (lj/B) = 0.287 0.687385 (max) 
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TABLE II. The maximal rotation velocities at R for differ- 
ent initial values of cr(0) for A = 10; cf. Fig. ^. 
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FIG. 1. The mass M ( — ) in units of (m P1 /m) and 

the particle number TV ( )in units of (mj^/m 2 ) as 

function of the central value (j(0) for different values of 



A := Xm^/Anm 2 = -5, 0, 5, 10 



FIG. 2. The redshift z depending on the self-interaction 
constant A calculated at the first maximum in the (M, a(0)) 
curve (see Figure |lj) which gives the last stable boson star 
configuration. 



FIG. 3. The redshift function 
self-interactions. 



for different 



FIG. 4. The rotation curves for the 
'normal' non-interacting boson star and with a large value 
of the constant A. Both curves are taken for the last stable 
boson star configuration with the largest mass value. 



FIG. 5. Rotation curves for different initial values 
<t(0) = 0.225,0.15,0.05,0.001 of A = 10; largest initial value 
has largest maximal velocity. 
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